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Abstract 

In the framework of a two-moment photo-hydrodynamic modehing of radiation transport, we 
introduce a concept for the determination of effective radiation transport coefficients based on the 
minimization of the local entropy production rate of radiation and matter. The method provides 
the nonequilibrium photon distribution from which the effective absorption coefficients and the 
variable Eddington factor (VEF) can be calculated. The photon distribution depends on the 
frequency dependence of the absorption coefficient, in contrast to the distribution obtained by 
methods based on entropy maximization. The calculated mean absorption coefficients are not only 
correct in the limit of optically thick and thin media, but even provide a reasonable interpolation 
in the cross-over regime between these limits, notably without introducing any fit parameter. The 
method is illustrated and discussed for grey matter and for a simple example of non-grey matter 
with a two-band absorption spectrum. The method is also briefly compared with the maximum 
entropy concept. 
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I. INTRODUCTION 



Excessive effort is required for modelling and simulation of radiation heat transfer in 
media with complex optical absorption spectra, like hot gases or plasma [l| . If the radiation 
model is part of a larger model for hot, compressible mixtures of various chemically reacting 
species consisting of complex ions, electrons, neutral molecules etc., and subject to transonic 
and turbulent flow, an exact treatment is not nearly possible. Applications range from arc 
physics in welding or electrical switching j^, 3, 4], atomic explosions, up to astrophysics j^. 
A frequently used approximate radiation model is photo-hydrodynamics, which is based on 
a two-moment expansion of the radiative transfer equation and a variable Eddington factor 

n 

(VEF) closure m. This concept can even be realized in a multiband framework, where the 

Q 

relevant quantities are decomposed according to their spectral properties [7|. 

The main problem of photo-hydrodynamic models is the optimal choice of the effective 
transport parameters, i.e., the effective absorption coefficients and the Eddington factor, 
which generally depend on the hydrodynamic and thermodynamic variables. The most 
prominent examples are the Planck mean absorption coefficient for optically thin media, the 
corresponding Rosseland mean fo. opfcaUy tMc. .nedia g a.d t.. constant EdCngton 
factor of 1/3 for isotropic radiation. 

Besides these special limit cases, the optimal definition of effective transport parameters is 
not straightforward j^, 1^. Accurate treatment of the general case is particularly important 
when radiation in the cross-over range between optically thin and thick limit dominates the 
physical behavior of a system, or if a medium is simultaneously transparent and opaque for 
different relevant wavelength bands. A simple approximate approach to solve that problem 
can consist in the construction of fitting expressions that interpolate between the limit 
cases. This has been done, for instance, by Sampson [11] for the absorption coefficient 
and by Kershaw 1^ for the variable Eddington factor. Another suggestion by Patch |l3| 
generalizes heuristically averages, which are exact for special cases. A very common approach 
is based on entropy maximization ^, 1^, [l^ . However, Struchtrup [l6| has pointed out a 
weakness of this procedure due to the neglect of the of radiation-matter coupling for the 
determination of the nonequilibrium photon distribution function, which plays indeed a 
main role for equilibration [l^]. One of the consequences is that already the Rosseland 
mean aborption coefficient in the near-equilibrium case is not correctly reproduced by a 
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two-moment photo-hydrodynamics with the maximum entropy closure. In this paper, we 
propose to use an entropy production rate principle, and we will show that the limit cases 
are correctly reproduced and the cross-over between them is provided in a natural way, i.e., 
without further model parameters. 



Maximization and minimization of the entropy production rate S have turned out to be 
powerful approaches for modelling many complex nonequilibrium systems [18||. We men- 
tion that whether the optimum of S* is a maximum or a minimum depends on the type 
of constraints 1^ and emphasize that entropy production optimization is in general not 
an exact physical law except near equilibrium, i.e., in linear deviation from equilibrium. 
However, it often provides useful approximate results even far away from equilibrium, pro- 
vided a predominance of strongly irreversible equilibration processes. It is also important 
for the discussion below that the method is not restricted to systems in (partial) local equi- 
librium, i.e., where the notion of (probably several) local temperatures can be introduced 
e.g., electron and ion temperatures in a plasma, light pencil temperatures for radiation 
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2l| . etc.). Entropy production optimization has been shown to be applicable to local 



nonequilibrium systems, for instance to Knudsen layers in material ablation or evaporation 
processes 
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23| . In such cases, the notion of entropy can still be defined [24 1. 



Entropy production of radiation has been discussed by Oxenius [20| and KroU 21j. Var- 
ious results related to entropy production principles in radiation have been reported. Essex 
25( 1 has shown that the entropy production rate is minimum in a grey atmosphere in local ra- 
diative equilibrium. Also this author has pointed out that a consideration of the interaction 
between radiation and matter is crucial because it contains the appropriate equilibration. 



27| discussed entropy 



i.e. entropy production, mechanism. Later on, Wiirfel and Ruppel 
production rate maximization by introducing an effective chemical potential of the photons, 
related to their interaction with matter. Finally, we mention Santillan et al. 28|| who showed 
that for a constraint of fixed radiation power, the black bodies are those which maximize 
the entropy production rate. 



This paper is organized as follows. In order to fix notation and to introduce the relevant 
quantities, in Sect. H the photo-hydrodynannc equations „e recalled. We will consider a 

system that is characterized by similar assumptions as in [13]. Scattering as well as photon 
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time of flight effects are assumed to be negligible, the matter is non-relativistic and in 
thermal equilibrium, and the ordinary index of refraction is unity. Section IIIII provides the 
basic result, i.e., Eq. (125]) determining the non-equilibrium photon distribution function. 
From this distribution function the transport coefficients can be calculated. In Sect. IIVI 
we show that the approach provides the Rosseland mean in the corresponding limit case. 
In Sect. |V] cases far from equilibrium are discussed. In particular, the emission limit, grey 
matter, and a simple artificial but illustrative example for non-grey matter are investigated. 
Our results are compared with the method of entropy maximization in Sect. IVI[ 



II. PHOTO-HYDRODYNAMICS 

We start from the Boltzmann equation (which is equivalent to the radiative transfer 
equation [1]) for the photon distribution function f{x,k) at location x and wave number 
k = ks, with direction vector s, | s |= 1. The function / is related to the spectral radiation 
intensity (radiance) / by / = {u^/c^)hf, where c is the velocity of light and u = ck is the 
angular frequency. Neglecting scattering, the Boltzmann equation for / reads Q (here, 
l|, K{k) includes the material density) 



in contrast to Ref. 



dJ + CS-Vf = -CK{k){f-f(^''^) . (1) 



Assuming that matter is in thermal equilibrium (at temperature T), the emission term 
contains the well-known Planck distribution 

^^''^^^) = ,1 1 ' (2) 

exp[nck/ ksi j — 1 

with the photon density of states y = 2/{2tt)^ including two polarization states. (For a 
generalization to nonequilibrium matter, e.g., a plasma with separate electron and ion tem- 
peratures, the emission term has to be replaced by the appropriate emission source function, 
and the entropy production of matter has to be treated in an appropriately generalized way.) 

In Eq. scattering is neglected, absorption and spontaneous as well as induced emis- 
sion are taken into account in the spectral absorption coefficient K,{k) = /t('^)(A;)[l — 
exp{—hck/kBT)]. 

It is well-known that in local equilibrium particle gas dynamics, from the Boltzmann 
transport equation hydrodynamic balance equations for mass, momentum, and energy bal- 
ance can be derived. For photons one can proceed in a similar way but with two serious 



restrictions. First, because the photon number is not conserved, a continuity equation anal- 
ogous to mass balance will not appear. The first moment is related to the energy, which 
is proportional to wave number for massless particles. Secondly, because the photon gas is 
generally not near equilibrium, one should in principle consider a large number of moments 



of general order and degree ly, |29|. Although we will restrict the number of moments to 



two, the entropy production method to be used is principally not limited to a specific num- 



ber of moments (cf. also SOj). 

In analogy to the P-N model 8|, an expansion can be defined in terms of spherical 
harmonics, 

fim = J d'kfik,e,^)kY;je,^), (3) 

where G and if are the zenith and azimuthal angle of the vector k, and the asterisk indicates 
complex conjugation. Yj^ are the usual spherical harmonic functions with indices / and m. 
Here, / gives the order of the moment. Due to the factor s in Eq. ([T]), the equation of motion 
for moments of order / is linked to moments of order / + 1. We introduce this notation, 
because for the analytical calculations below, it is sometimes convenient to use spherical 
harmonics. 

In Cartesian coordinates, the first three moments fim are associated with energy (/ = 0), 
momentum (/ = 1), and radiation pressure tensor (/ = 2). Energy (photon energy hu = hck) 
and momentum (photon momentum hk) densities are defined by 

e(f) = [ £kfhck , (4) 



p'(f ) = J d^k fhsk . (5) 
The relations of energy and momentum to the moments defined in ([3]) are, for instance 



e = hcy/AnfQQ and = h^jAn/Sfio, respectively. Multiplication of Eq. ([T]) with tick and 
hks and integration over momentum space gives 

dte + c^V-p = Pe , (6) 
dtp + V ■ n = 4 , (7) 

where 11 is the pressure tensor of the photon gas, having the components 

n-mn — / d^k SmSnffick . (8) 



The radiation pressure is S^^]^n„„/3 = e/3. The source terms of Eqs. ([6]) and ([7]), which 
mediate the interaction of radiation with matter, are 

Pe = -he' j £kkK{k)U - f^'''^) , (9) 

Pp = -he I d^kkK{k)f . (10) 



The radiative heat production w{x), which will appear as a source term in the energy balance 
equation for the matter is given by 

W = -Pe . (11) 

For a closure of two-moment expansion, the pressure tensor 11 has to be expressed as a 
function of energy e and momentum p. Based on tensor symmetry arguments and on the 



trace requirement, one can show that 11 has the general form pj| 

Ilmn = e I ^ dmn H ^ ^ I , (12) 

where x(e?p) (with p =| p \) is the so-called variable Eddington factor (VEF) describing 
radiation anisotropy. If radiation is isotropic (e.g., in the purely diffusive limit) 11 must 
be proportional to the identity matrix, which implies x = 1/3- For a beam (e.g., p in z- 
direction), however, x = 1 must hold. Therefore, x allows to describe non-diffusive, directed 
radiation. The pressure tensor is related to the second order moments f2m- For instance, 
assuming p in z-direction, Eq. ([8]) and I20 = \/5/167r(3 cos^(0) — 1) imply the relation 



nss = e X = ^ + y ^/ic/ao j , (13) 

which will be needed below. Effective absorption coefficients be formally defined by 

expressing the source terms and (ITUl) by 

Pe = c/€e(e(^'')-e) , (14) 
Pp = -eKpp . (15) 

The effective absorption coefficients Ke,p{e,p,T) generally depend on e, p, and T. Here, e^^'^'' 
is the equilibrium energy density 

c 

where a = 27r^fcg/15c^/i^ ~ 5.67 10^^ W/m^K^ is the Stefan-Boltzmann constant. The 
energy radiated by a blackbody per time and area is given 



Equations and ([7]) define the photo-hydrodynamic equations for tlie variables e and p 
of the photon gas. In order to solve a complete radiation problem, Eqs. ([6]) and ([7]) must be 
solved together with the hydrodynamic equations for the matter. This work focuses on the 
problem of the determination of the yet unknown effective transport coefficients Kg p and x- 
For this, the distribution f{x,k) must be known. Often simple approximate averages are 
used involving the equilibrium function f^'^'^\ like the above mentioned Planck mean, 

he I'd^kKkf^^'^^ 
'^^ = eM ' 

and Rosseland mean, 

_J^klfd^J^ 

^° jd^kK-^k^dkf^-D ■ ^ ' 

While the former is a weighted spectral average of the absorption coefficient n{k) that 
depends on the wave number k, the latter is the inverse of a weighted spectral average of 
the inverse absorption coefficient The two cases strongly differ as the Planck limit is 
dominated by wave number bands with large n{k) values, while the Rosseland limit is dom- 
inated by small K.{k). In the associated limit cases radiation is isotropic and the Eddington 
factor, being the ratio of pressure and energy density, equals 1/3. 

Prior to the calculation of / and the effective transport coefficients with entropy produc- 
tion minimization, we mention that the issue of photo-hydrodynamic boundary conditions 
or e and p will not be discussed in this paper. We refer the reader to appropriate literature 
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III. MINIMIZATION OF THE ENTROPY PRODUCTION RATE 

In the following, we determine /, and from this xi^iP) ^"^^ i^e,p{G,p), by minimizing the 
total entropy production rate S under constraints of fixed values of e and p (we always 
assume p in z-direction). First, we derive an expression for the^elevant part of the entropy 
production rate. The radiation entropy is given by (see, e.g., 

SUx) = -kB j d'k (^f \n{f/y) - {y + f) ln(l + f/y)j . (19) 

The part of the entropy production rate of the photon gas alone is obtained by differentiation 
of Eq. (HM with respect to time and subsequent replacement of dtf with the help of Eq. ([T]). 
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This leads to an equation S^S'rad + Vg^^^^ = S, with qs^^^ being the entropy flow density, and 
the local entropy production rate 

S(f ) =kBj d'k In {y/f + 1) c/t(/(^'') - /) . (20) 

The total matter-radiation system includes the hydrodynamic equations for the matter. 
The entropy production rate of matter contains a matter-specific radiation-independent 
part, which is constant under variation of / and is thus not of interest here. Additionally, 
the energy balance of matter contains the power density Eq. ( ITTi) . which can be associated 
with a local entropy production rate, w/T, where w is obtained from Eq. ([9]). The /- 
dependent part of the total local entropy production rate of the radiation-matter system is 
thus S = T, + w/T. Using Eq. ([2]) for replacement of the temperature T leads finally to (see, 
e.g., M) 

Six) = -ksc I d^kln [ yffitll ) <f - • (21) 

Optimization of S subject to the constraints of fixed p and e implies that the variation of 

e(f ) - j <fk fhck^ - ■ (^p{x) - j d^k fhk^ (22) 

with respect to / vanishes. It follows from the form of Eq. fl2T]) that the optimum is 
a minimum, which we confirmed also numerically. In Eq. (122|) dimensionless Lagrange 
parameters Ae and \p have been introduced. Variation gives 

/-/'"' . ( y/f + l \ . k r- k 




This is the central result of this paper. The appearance of K.{k) expresses radiation-matter 
interaction as the entropy generation process. 

The radiation heat transfer model proposed in this paper can now be summarized. Equa- 
tion (l23ll delivers implicitly /(/c, T, Ag, Ap). Using this /, from the constraints (jlj) and ([5]) 
one can derive e(Ae, Ap,T) and p(Ae, Ap,T), as well as fi;e(Ae, Ap,T) and Kp(Ae, Ap,T), all of 
them still dependent on Ae,p. These Lagrange parameters are functions of e,p and T by 
virtue of e = e(Ae,Ap,T) and "p = p{Xe, \p,T). Appropriate replacement leads then to the 
wanted functions x(e,p, T) and Ke^p{e,p,T). 

Although this procedure can be cumbersome, for a given material (gas, plasma, etc.) 
it can be done once, and the results can then be stored in look-up tables or described by 
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appropriate fit functions, which can be used further in the simulation of the hydrodynamic 
approximation. In the following section we discuss specific cases that are simple enough for 
an analytical treatment. 

Before we discuss the results, we anticipate a remark on the maximum entropy approach 
(see e.g. jo], Q), where S'rad takes then the part of S in Eq. (12^ . which has to be optimized. 
Because, in contrast to S in Eq. (12T|) the entropy S'rad in Eq. (fT9!) does not depend on the 
absorption spectrum K{k), the resulting ME distribution function /me does neither. Hence, 
in contrast to our approach, the ME approach does not take into account radiation-matter 
interaction explicitly at the level of the photon distribution function. 

As a side remark, we mention that irreversibility, due to radiation-matter coupling, again 
enters in both approaches at the hydrodynamic level, i.e., when Eqs. (E]) and ([7]) are solved 
together with hydrodynamic equations for the matter. Therefore, eventually both methods 
exhibit irreversibility. But ME assumes that radiation is in a conditional maximum entropy 
state, while our approach goes one step further by considering nonequilibrium states away 
from the entropy maximum. 



IV. RADIATION NEAR THERMODYNAMIC EQUILIBRIUM: ROSSELAND 
LIMIT 

If the deviation from equilibrium can be approximated by linearization ('near equilibrium' 
or 'weak nonequilibrium' of the radiation), one can expand Eq. fl23|) in a Taylor series with 
respect to 5/ = / — f^^"^ and keep only the leading order terms. In this section we will use 
spherical moments ([3]), in order to show that the approach yields the Rosseland limit up to 
arbitrary order / of moments. Furthermore, we make an expansion to second order in 6f, 
which yields for the VEF the first nontrivial order in p beyond the constant value of 1/3. 

Vanishing variation of entropy production with constraints of fixed moments fim reads 
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S' - ^ kscXim ifim- / d'^k kfYi*^ 



. (24) 



Here, 5 denotes the variation with respect to / and Xim-, the constant factor /c^c in the con- 
straints is used to obtain dimensionless Lagrange- multipliers Xim with indices I and m, and 
the sum extends over all moments that are to be constrained. By expanding the expression 
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in 5f up to second order and solving in second order of the Lagrange multipliers, we obtain 

= -f''\y + f'"^) /i + 3(2/ + 2/(^'^))/(^'^)(y + ^ (25) 

with ji = {k/2yK{k)) ^ A^m^m) which must be real. The unknown Lagrange multipliers Xim 
must be determined from the constraints = / cPk k{f^^''^ + 5f )Yi^. Insertion of 5f from 
Eq. (j25l) gives 

flm = SmoSlofoQ'^^ + CiA;*^ ~^ y ^^^^l*rn^ ' (^^) 

with /(So"^ = e^'^^V^cVivr and 

We used f^^''\y + f^^''^)hc = —ykBTdkf^'^'^\ which simplifies both terms on the right hand 
side of Eq. fl25|) . Solving Eq. fl26l) for the Lagrange multiplier A;*^ and inserting into Eq. 
( !25|) . gives to leading order in A 

.._kBTd^ k f^^^^^ 

- 2hc W) ^ C[ • ^^^^ 

Generalizing f[T^ . f[T^ . effective absorption coefficients for mode I and m can be defined by 

Plm = —l^lm [flm — ^mO^lo) ■ (30) 

Inserting (l29l) into the source term Pim = —he f d^kknf [see Eqs. ([9]), (fTOl) ]. we obtain 

_ Jldkk^d^f^ _ 

In this approximation, the absorption coefficients do not depend on the order of the moment 
and are equal to the Rosseland mean f|T8|) . Hence, equilibrium. This 

is not very astonishing, because the Rosseland average is the appropriate effective radiation 
absorption coefficient near local thermal equilibrium and, in this limit, the entropy produc- 
tion rate principle holds exactly. 
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Let us now calculate the Eddington factor x near equilibrium. Only constraints with 
/ < 1 have to be used. Aligning the 2;— axis of the coordinate system along the momentum 
p, the Eddington factor is given by x = nas/e, and all moments with m 7^ vanish. We will 
thus ignore the index m in the following (e.g., /o = /oo, /i = /lO; etc.). The two equations 
for Z = 0, 1 can be solved for A; up to second order in deviation from equilibrium: 



Ao 
Ai 



/i 



Co 



c 
1 



'0 

2Cr 



(/o-/o^^^^ + /i)V(/o-/o^^^^-/i 



{fo-ft^ + fiy-{fo-ft^-fi 



(32) 
(33) 



This allows to calculate from Eq. fl25|l up to second order 6f and /20 needed in Eq. f|T3l) . 
One finds /20 = (^2/1^ /"s/SvrCf . Using = p = CA/47r/3/i, the Eddington factor becomes 
up to quadratic order away from equilibrium 



v ^ — 
X 3 



he 



Co 



75Ct 



cp y 



(34) 



Thus, near equilibrium, the Eddington factor is 1/3 as it must be, and the VEF grows 
quadratically with p. The curvature depends on K,{k) unless the medium is grey, as will be 
discussed in Sect. IV B[ 



V. GENERAL NONEQUILIBRIUM RADIATION 
A. The emission limit 

If the photon density is so small that absorption can be neglected, the emission approx- 
imation can be applied. Formally, this refers to the limit / <^ f^^'^\ or e/e*^^^-* 0, such 
that the source term on the right hand side of Eq. ([T]) becomes CK{k)f^^'^\ It is clear from 
Eq. ([9]) that the matter is then cooled, with a cooling power density w = — CKpie*-^'^-' that 
contains the Planck mean (fT7|) for Kg- In order to determine Kp, we must go to first order in 
/ beyond the Planck limit. Expansion of Eq. fl23l) up to leading order of f / f^^"^^ gives 

/ = — ^^i^/(-)(fc) . (35) 

Substitution of / in Eqs. (jlj) and ([5]) leads to implicit equations for e(Ae, Ap) and p{\e, Ap). 
They can be solved in leading order of Ap corresponding to the case of small momentum 
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p or almost isotropic radiation. One finds Ae = hcR/e and Xp = —Shc^kp/e^ with k = 
J d^hnf^'^'^K Replacement of / in Eq. (fTOj) and using these results leads to Pp = —n^cp/k, 
with = j d^kK? f^^'^\ In summary: 



Ke = Kpi (36) 

Kp = ^ . (37) 

r\j 

Energy and momentum relaxation have different effective relaxation constants in the limit 
of small e for non-grey matter. From Eq. fl35l) one obtains with Eq. ([8]) the Eddington factor 
X = 1/3 in this limit, i.e., in the emission approximation radiation is also isotropic. 



B. Grey matter 

Grey matter is characterized by wave number independent spectral absorption, /t(A;) = 
const. The effective absorption coefficients and the nonequilibrium distri- 

bution / does not depend explicitly on the k value. This follows from the proportionality 
of the entropy production rate to /t in Eq. ( !23l) . However, note that in the framework of 
the hydrodynamic equations ([6]) and (I?!), / depends implicitly on the n value via its e and 
p dependence. 



We have numerically calculated / for grey matter. In Fig. [T] a plot of S,^f is shown as a 
function of ^ = hck/ksT for different values of e at p = 0. One observes that it is mainly 
the photon occupation number, which increases for higher energy e, while the shift of the 
wave number at maximum / towards larger values is comparatively weak. 

As the momentum p increases, however the wave number is expected to be more affected. 
Figure [2] shows the distribution as a function of wave number for different values of the 
direction cosine. In this example, pc = e^'^'^^ and e = 2e^^'^\ which roughly means that half 
of the energy is thermal and half is "ballistic". The function C,^f at cos(O) = is thus 
very close to the equilibrium distribution for e = e^'^'^\ which is also plotted in the figure. 
Besides a change in their number density, the photons parallel (cos(B) = 1) and anti-parallel 
(cos(O) = —1) to the beam are shifted towards higher and lower wave numbers, respectively. 

The variable Eddington factor (VEF) obtained from our results turns out to be close to 
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the heuristic proposal by Kershaw [12], Xk = (1 + 2(cj9/e)^)/3. In particular, as required 
X = 1/3 for p = and x 1 for PC e. Figure [3] shows a plot of x as a function 
of cp/e for various cases, xk seems to be a reasonable approximation for the VEF. For 
grey matter the integrals (I27p and fl25]) can be calculated analytically by partial integration. 
Equation (1341) implies then for the VEF near equilibrium (i.e., e ~ e'-^'^^ and pc <^ e) 
x(e,p)^(l + (3cp/2e)2)/3. 

The results obtained from the entropy production method are not only exact in the 
weak nonequilibrium limit, but provide the correct VEF in the limit of directed radiation 
{p e/c). This supports the conjecture that the entropy production method, in the present 
context, can serve as a useful approximation even far from equilibrium. 

C. Non-grey matter 

As a simple example with mainly illustrative purpose, we consider now an artificial two- 
band absorption coefficient. We will use the non-dimensional quantity ^ instead of the wave 
number k and consider a step-function absorption-spectrum of the form k{S,) = 2ki for ^ < 
and k{^) = Ki for ^ > C,c, where ni is constant (cf. Fig. H]). This artificial spectrum k(^) 
describes matter-radiation interaction that is stronger for long than for short wavelengths. 
To be concrete, = 4, which gives Kpi = 1.6 ki, kro = 1-26 ki, and n'^/k = 1.89 ki. 

One might have expected that for e — > (emission limit) the value of Kp becomes equal 
to K,{C, ^ 0) = 2/ti, which is obviously not the case. The value Kp = 1.89 ki in this limit can 
be understood with the help of Eq. (!35l) . which yields / ~ f^'^'^^K(k)e/}ihck. Hence / is a 
given function of k (or and is proportional to e. This implies that Kp is a fixed weighting 
of the spectral absorption coefficient. 

The function ^^f obtained from the minimum entropy production principle is shown in 
Fig. mfor p = and different e/e*^^''^ values. A similar qualitative dependence of the photon 
occupation number on e/e'-^'^^ appears as for grey matter. Furthermore, the step in k{^) 
leads to a step in /, such that the nonequilibrium distributions are closer to equilibrium for 
C, < ^c- This pull of the long-wavelength part towards equilibrium can be understood from 
the larger k{^) value in this region. Indeed, stronger interaction of radiation with matter 
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leads to stronger equilibration. 



The VEF for this special case of a stepwise absorption spectrum is shown in Fig. O The 
VEF has increased as compared to the grey medium, hence the deviations from Kershaw's 
VEF are a little bit larger. The increase of x at fixed cp/e might be understood qualitatively 
as follows, if one interprets x as a measure for the radiation pressure in direction of momen- 
tum p. Because the considered absorption spectrum leads to enhanced equilibration of the 
long wavelength photons, there must be an increased amount of short wavelength photons 
contributing to the momentum p, which leads to the larger x 

The normalized effective absorption coefficient Kg/ ki is shown in Fig. [HI as a function 
of e/e^^''-' at p = 0. In the limits e and e = e^^'^\ the Planck and Rosseland mean 
are obtained as one expects. For e ^ oo it holds Ke{e,p = 0) — ki, as most photons 
will populate the short wavelength band. For finite p, it is inconvenient to discuss Ke{e,p) 
because the zero of Pe{e,p) is no longer given by e = e'-'^^^ but shifts as a function of p. (Note 
that finite p is always associated with nonequilibrium.) An effective Kg, however, can still 
be defined as in Eq. (Q, but has a pole at e*-^''-' and a zero shifted away from e'-^'^^ This fact 
is illustrated in Fig. [TJ which shows that Pe{e,p) starts to deviate from the value at p = 
as p increases (particularly the shift of P^(e^'^'^\p) from zero). Although equiva^ 
be more convenient to work with or with representations as in, e.g., Ref. 
Kg, if general simulations of a fully coupled hydrodynamic radiation problem are performed. 



ent, it might 
3] than with 



Figure [8] shows the effective absorption coefficient K,p{e,p)/ ki as a function of pc/e for 
different e/e^^'^-'-values. Again, for small p, Kp is large at small e (emission approximation, 
where it goes to the value 1.89), while near equilibrium the Rosseland mean is reproduced. 
At large e, Kp goes to /ti because the distribution function / extends to higher and higher 
wave numbers. The same behavior occurs when cp/e becomes large, i.e., Kp approaches ki, 
as can be concluded from the dotted curve in Fig. [21 which indicates that the distribution 
function of directed photons involves large wave numbers. 



14 



VI. COMPARISON WITH THE MAXIMUM ENTROPY METHOD 



In this section we compare our method with the well-known ME approach, where the 
radiation entropy is maximized. Especially, effective absorption coefficients and VEF are 
compared for the different explicit cases treated in the previous sections. The basic drawback 
of ME is the explicit independence of / on K,{k). The implicit dependence on k, via radiation- 
matter coupling at the hydrodynamic level (i.e., via an e and p-dependence due to the ME 
constraints) is not able to spectrally resolve the absorption properties in the distribution 



function. Indeed, one can easily calculate /, which is given by 
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/me — ~ — ~ (38) 

exp {hck{Xe — Xps)) — 1 

with e and p dependent Lagrange parameters. At this level, information about the absorp- 
tion spectrum is not contained. 

On the other hand, according to Eq. (l23l) . the distribution function / obtained from 
entropy production minimization does explicitely depend on K,{k). In subsection IV C[ it 



has been shown that there is a significant and easily understandable effect of wave number 
dependent K{k) on the distribution function, and thus also on the effective absorption coef- 
ficients and on the VEF. 

First, we note that because of the /t-independece of /me, the associated VEF is a unique 
function of e and p, 

5 4 f 3^ 

which is the well-known Levermore-Eddington factor [16]. It is also plotted in Figs. [3] and 
|5] for comparison. For grey matter the Eddington factors obtained by our approach and 
by ME are different, but probably similar enough to justify the use of Xme- For non-grey 
matter the deviations increase. 

The difference in the effective absorption coefficients Ke^p is more important. In the limit 
case of weak nonequilibrium, where entropy production minimization reproduces the correct 
Rosseland mean, ME gives (i.e., up to leading order in /me — f^^''^) 

_ j d^kKk'^dkf^^"^ 
- -fd^kk^dji^ ' 



which is clearly different from the correct Rosseland mean (fT8|) (cf. 16j). Note that the two 



methods yield different results for the transport quantities in equilibrium limit, because the 
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deviation from equilibrium, Sf = f — f^^''\ is different. Although the zero'th order term 
is the same and 5f goes to zero, it is 5f which determines the transport coefficients. 
As an example, we have added the curve for kme in Fig. [61 While the Planck limit (e — > 0) 
is correctly obtained by ME, there is a difference of about 10% to the Rosseland limit 
(e = e^'^) for our simple example. As has been shown in Ref. 16|, the ME approach is 
able to reproduce the Rosseland mean only if one considers a larger number of moments, 
which goes beyond a two-moment photo-hydrodynamic description of the photon gas with 
variables e and p. 



A remark on the K-dependence of our approach is in order. One might conclude from 
Eq. fl2T]) that in the case of vanishing n the minimum entropy production rate approach 
is not defined, while the maximum entropy approach is defined because n does not appear. 
However, zero k means absence of equilibration. But maximum entropy requires that equi- 
libration to the associated maximum entropy state already occurred. Because for zero k 
there is no process that drives the photon gas towards the associated maximum entropy 
state, both approaches encounter a similar conceptual problem. 



VII. CONCLUSION 



The description of radiative heat transfer with photo-hydrodynamic equations requires 
a closure procedure to restrict the number of coupled equations. If radiative energy den- 
sity e and momentum density p are to be considered, effective absorption coefficients and 
a variable Eddington factor (VEF) have to be determined as a function of e and p. Using 
the minimization of the (local) entropy production rate of radiation and matter, a closure 
procedure is introduced that yields the photon distribution function, which depends on the 
spectral absorption coefficient, and which allows for the calculation of transport coefficients 
in general. 

It has been shown that the derived expressions are exact near equilibrium, in the emission 
limit, and correctly describe directed radiation with a VEF x = 1- Whereas the first fact is 
not very surprising as minimum entropy production is known to be exact near equilibrium, 
the other correct limiting behavior demonstrates that entropy production optimization can 
provide sound and useful results also far from equilibrium. 
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For the general case, effective absorption coefficients and a variable Eddington factor are 
found, which reasonably interpolate between the limiting cases, notably without introduc- 
tion of any fit parameter. The variable Eddington factor for grey matter is close to the one 
proposed by Kershaw [l^, and not far from the Levermore-Eddington factor, which can be 
derived from a maximum entropy approach. 

In contrast to our approach, the distribution function obtained from the maximum en- 
tropy (ME) method does not contain the absorption spectrum. Hence, ME cannot describe 
frequency resolved absorption effects on the distribution function. An example of such an 
effect, which is taken into account by our method, is the enhanced photon equilibration in 
wave number bands with larger absorption as discussed in Sect. IV CI Another difference 
is that in the two-moment photo-hydrodynamic framework, ME is unable to reproduce the 
Rosseland mean absorption coefficient in the equilibrium limit. 

Roughly speaking, the maximum entropy approach is a zero order approximation in the 
sense that radiation has equilibrated to a conditional maximum entropy state. Our approach 
can be understood as a first order approximation by assuming that the radiation is away 
from equilibrium and equilibrates according the minimum entropy production method. To 
summarize, we conjecture that this method may serve as a useful approach in various radi- 
ation heat transfer problems. 

A benchmarking and critical discussion of specific application examples will be necessary 
in the future. Furthermore, the feasibility of improvements in the framework of multiband 
extensions or by taking into account higher order moments should be investigated. 
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FIG. 1: Function /^^ for p = as a function of ^ = hck/kBT (in analogy to the spectral power 
density) obtained from the entropy maximization method applied to a grey body, e = e^^'^^ (solid, 
Planck function), e = 2e(^^) (dotted), and e = O.Se^^^^ (dashed). Larger energy implies a consider- 
ably higher photon occupation number and only a weak shift to larger wave numbers. 
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FIG. 2: Function /^^ as function of ^ = Hck/kBT and for different values of the cosine of the angle 
between k and p. The value cos(B) = 1 corresponds to photons parallel to p, while cos(0) = — 1 
corresponds to anti-parallel photons. Here, the energy is roughly equipartitioned between ballistic 
and diffusive motion, i.e., e = 2e^'^'^\ p = e^^'^^/c; symbols indicate the equilibrium distribution. 
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FIG. 3: VEF of a grey body for the different cases listed in the legend. Deviations from xk are 
below the percentage range, deviations from the Levermore-Eddington factor xme are larger. 
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FIG. 4: Top: Step-function absorption coefficient. Bottom: Function f^^ for different eje^^'^- 
values. The change in k at ^ = 4 imphes a change in the distribution function. Higher values 
pulls the distribution function / stronger towards the equilibrium distribution j^^'^ because of 
stronger radiation-matter interaction. 
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FIG. 5: Same as in Fig. 3, but now for the step-function absorption coefficient of Fig. 4. The 
deviations from xk are a few percent, while deviations from xme are twice as large. 
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FIG. 6: Effective absorption coefficient Kg (divided by ki) at p = for the step function absorption 
coefficient of Fig. 4. obviously takes the Planck mean value (dashed-double-dotted) at e = 
(emission limit) and the Rosseland mean (dashed-dotted) at e = e^^^^ (equilibrium radiation). For 
large energies, Kg will converge to ki, i.e. to the n{k) value at large wave numbers populated by 
most of the photons. Entropy maximization leads to a mean absorption kme (dotted) with the 
correct Planck mean at e = 0, but with a deviation from the Rosseland mean at e = e*^^^). 
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FIG. 7: Energy loss function as a function of cp/e for different eje^^'^ values. This figure 
illustrates that the zero of Pe(e) is located at e^^'^ only for equilibrium conditions, i.e., p = 0. For 
finite p the zero slightly shifts, i.e. Pe.ie^^'^^ is finite. 
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FIG. 8: Effective absorption coefficient Kp (divided by ki) as a function of cp/e for the step function 
absorption coefficient of Fig. 4. Tfie dashed-dotted and the dashed-double-dotted lines indicate 
the Rosseland and Planck means, respectively. 
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